Study of Doubly Heavy Baryon Spectrum via QCD 

Sum Rules 



Liang Tang a , Xu-Hao Yuan", Cong-Feng Qiao a ' & and Xue-Qian Li 

a) School of Physics, Nankai University, 300071, Tianjin, China 
b) Department of Physics, Graduate University, the Chinese Academy of Sciences 
YuQuan Road 19A, 100049, Beijing, China 



In this work, we calculate the mass spectrum of doubly 
heavy baryons with the diquark model in terms of the QCD 
sum rules. The interpolating currents are composed of a 
heavy diquark field and a light quark field. Contributions 
of the operators up to dimension six are taken into account 
in the operator product expansion. Within a reasonable 
error tolerance, our numerical results are compatible with 
other theoretical predictions. This indicates that the diquark 
picture reflects the reality and is applicable to the study of 
doubly heavy baryons. 
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1 Introduction 



The considerable success of quark model in interpreting a large amount of hadronic ob- 
servations has convinced people its undoubted validity for many years. In the quark model, 
hadrons are constructed according to two configurational schemes: mesons, consisting of 
a quark and an antiquark (qq); and baryons, consisting of three quarks (qqq). Right after 
the birth of the quark model, the diquark model was proposed where two quarks constitute 
a color-anti-triplet which behaves as an independent object in the baryon. In Gell-Mann's 
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original paper on the quark model, he discussed the possibility of the existence of free 
diquarks [I]. The concept of diquarks, has been established in the fundamental theory, 
and has been invoked to help illuminating a number of phenomena observed in experi- 
ments El lU El El E] ■ The systems composed of three quarks should be described by the 
Faddev equations, but since there are three coupled differential equations, solving them is 
extremely difficult. As a matter of fact, the three-body problem is still an unsolved sub- 
ject even in classical physics. It is tempted to consider the diquark-quark structure which 
turns the three-body system into a two-body one, and the three Faddev equations then 
reduce to single equation (no mater relativistic or non-relativistic). Thus the problem is 
greatly simplified and solution concerning baryon physics is obtained. However, for the 
baryons which are composed of three light quarks, the three Faddev equations have the 
same weight, so a problem emerges right away, namely which two quarks are combined to 
compose a diquark while the rest one moves independently. It seems to be an unbeatable 
difficulty. However, recently the topic on diquarks revives, for it may bring up some direct 
phenomenological consequences. Especially, when there are two heavy quarks in a baryon, 
they may constitute a relatively tight structure, a diquark. A diquark has the quantum 
numbers of a two-quark system. For the ground state, a diquark has positive parity and 
may be an axial- vector (5 = 1) or a scalar (5 = 0). According to the basic principle of 
QCD, for the two quarks residing in a color anti-triplet, the interaction between them is 
attractive. 

Baryons containing two heavy quarks are important and intriguing systems to study 
the quark-diquark structure of baryons. The two heavy quarks (b and c) can constitute 
a stable bound state of 3, namely, a diquark which serves as a source of static color field 
for light quarks [8]. The SELEX Collaboration reported the first observation of a dou- 
bly charmed baryon via the decay process — > A+K~tt + , by which its mass of 3519± 
lMeV/c 2 [H] was determined. Later, this baryon was confirmed by the SELEX Collabora- 
tion through the measurement of a different decay mode — > pD + K~ , whose final state 
involves a charmed meson JTU]. However, both the BABAR and Belle Collaborations did 
not observe this state in e + e~ annihilation experiments [HI [T2]. This may be due to the 
distinct beam structures of the two types of experiments, and the reason is worthy of 
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further and careful studies. 

In the theoretical aspect, there have been numerous works in studying the doubly 
heavy baryons [131 HH EE HHJ E] • All these works concern the dynamics which results in 
the substantial diquark structure. 

Although the QCD is proven to be an undisputablly valid theory about strong inter- 
action, the non-perturbative QCD which dominates the low energy physics phenomena 
has not been fully understood yet. Among the the theoretical methods in dealing with 
the non-perturbative effects, the framework of the QCD Sum Rules which is indeed a 
bridge between the short- distance and long-distance QCD as initiated by Shifman et 
al. [18], turns out to be a remarkably successful and powerful technique for computing 
the hadronic properties. Recently, a number of works have been worked out to inter- 
pret the newly observed mesonic resonances within the framework of the QCD sum rules 
[T9l [20j [21] . Meanwhile with the QCD sum rules, a few works were performed in studying 
the mass spectrum of doubly heavy baryons [221 [231 12U [251 [26]. In those studies, the 
authors calculated the correlation function of baryonic currents composed of quark fields 
by virtue of the operator product expansion (OPE). 

Since the correlation of the two heavy quarks is strong, they are tempted to be bound 
into a diquark which can be regarded to manifest independent degrees of freedom in the 
baryon. In this work, we no longer treat the two heavy quarks as independent constituents, 
but a combined sub-system— diquark which behaves as a component of doubly heavy 
baryons, and the corresponding field is denoted by a new bosonic symbol $ with a mass 
mjj. In fact, this picture was recently proposed in Ref. |27j. Then, in this tentative model 
for calculating the mass spectrum of doubly heavy baryon systems, the light quark q 
(q=u,d,s) orbits the heavy diquark which is a tightly bound QQ' (Q=c,b) pair. The 
application of the diquark can simplify the interpolating currents which are important for 
obtaining the baryon spectrum in the QCD sum rules. The spin-parity quantum number 
of a ground-state diquark is either + or 1 + . The former, along with a light q, can form 
the state with J p = | + ; the latter can form not only the state with J p = | + , but also 
J p = | + . That is to say, using the model of the diquark and the QCD sum rules, we can 
study the doubly heavy baryons with spin-parity J p = | + and J p = | + . 
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The content of the paper is arranged as follows. In Sec. II we derive the formulas of 
the correlation function of the interpolating currents with proper quantum numbers in 
terms of the QCD sum rules. In Sec. Ill, our numerical results and relevant figures are 
presented. Section IV is devoted to a summary and concluding remarks. 

2 Formalism 

The method of the QCD Sum Rules is starting with choosing proper correlation func- 
tion (or Green's function) both at the quark-gluon level and the hadron level. The corre- 
lation function for the doubly heavy baryons reads 

mn-ij^wvmow). (i) 

Considering the spinor structures of baryons, the correlation function has the Lorentz 
covariant expression as follows [221 El] : 

n(g 2 ) = ^n 1 (g 2 ) + n 2 ( g 2 ). (2) 

For each invariant function of Ui(q 2 ) and n^g 2 ) in the doubly heavy baryons one can 
obtain a sum rule. 

Following Refs. [23l [24] and based on the diquark model, the interpolating currents, 



which play a crucial role in our analysis, are chosen to be 

J(x) = $*(x) T£ q a {x), for spin 1/2 baryon; (3) 

J*(x) = ^{x)T k q a (x), for spin 3/2 baryon; (4) 

J'(x) = & a (x) T fc q a (x), for spin 1/2 baryon, (5) 



where $t(x) and Q a (x) are axial vector and scalar diquarks, respectively. The interpo- 
lating current J corresponds to Sqq/ and Qqq>, corresponds to Sqq, and Qqq' anc ^ J' 
corresponds to Hqq/ and Q'qqi respectively with Q,Q' = c, b. The concrete definition of 
r£ and Tfc are presented in Table [TJ 

On the phenomenological side, the correlation function is expressed as a dispersion 



4 



Table 1: The choice of and IV The index D in J D D means the diquark. ${qq'} denotes 
the axial vector diquark and $[qq'] denotes the scalar diquark respectively. 
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"QQ' 
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1 + 
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0+ 
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integral over a physical regime, 

^ + Afg f°° p h (s) 

where Mjj is the mass of the doubly heavy baryon, A# is the baryon coupling constant 
and p h (s) is the physical spectral function of the continuum states. When we attain the 
above expression, the summing relations for the Dirac and Rarita-Schwinger spinors have 
been used, namely, for spin-3/2 baryons the numerator of the first term in Eq.([6]) should 
be replaced by the proper Lorentz structure 



(4 + M H )(g^-- W+ ^ 3M|; . 

With the operator product expansion (OPE), the correlation function Tli(q 2 ) (i=l or 
2) can be written as: 

6 

u t ( q 2 ) = nfV) + n r d ' dim (q 2 ) • (7) 

dim=3 

Here, "pert" , "cond" and "dim" refer to perturbative QCD calculation, the quark or gluon 
condensates, and the relevant condensate dimensions, respectively. nf ert (g 2 ) is obtained 
by taking the absorptive part of the Feynman diagram A, and n^ ond,dim (g 2 ) represents the 
contributions from various condensates. In this work, we consider the condensates up to 
dimension six, as people usually do in the literature. 
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Figure 1: Feynman diagrams. 



The Feynman diagrams contributing to the correlation function of the doubly heavy 
baryons are displayed in FigJU and the gluon-diquark vertices are shown in FigfSJ In 
Ref . |28] . the effective gluon-diquark vertices were given, which will also be used in our 
later calculations, so we just copy them below: 



SgS = ig s t a ( Pl +p 2 ),F s (Q 2 ), 

VgV = -g s t a {g a p(p 1 +p 2 )i 1 - g m [(l + ^px- K V p 2 ]p 
- g^[(l + k v )p2 - KvPi] a }F v (Q 2 ) ■ 



(8) 



(9) 



Here, Q 2 =—(pi —P2) 2 , g s = y/^ot s denotes the QCD coupling constant, k v is the anoma- 
lous (chromo) magnetic moment of the vector diquark and t a (= X a /2) is the Gell-Mann 
color matrix. Furthermore, Fs(Q 2 ) and Fy{Q 2 ) are the diquark form factors. 

The scalar diquark's propagator is p2 _^i +ie , and the axial- vector diquark's propagator 



is 



* 9v 



p^—mi+ie 



. Since the diquarks are made up of two heavy quarks, according to the 



general rule, their condensates are negligible [29J. 

The Feynman diagrams are computed not only with the regular QCD Feynman rules, 
but also with the effective vertices displayed above for point-like diquarks. For taking 
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Figure 2: Gluon-Diquark vertices. The wave line represents the gluon, the dashed repre- 
sents the scalar diquark, and the black line represents the axial vector diquark. 



into account the composite nature of diquarks, phenomeno logical form factors {Fg(Q 2 ), 
Fv(Q 2 ) ) are introduced. The authors of Ref. [28j suggested the vertex functions as fol- 
lowing: 



Fs(Q 2 ) 

5s 
5 S 

Fv(Q 2 ) 

5 V 
5v 



6 s n2 



Q 2 s 



Qs + Q 2 

a s (Q 2 )/a s (Q 2 s ) if (Q 2 > Q 2 S ) 
1 if (Q 2 < Ql) ; 



(10) 



Qv 



Q 2 v + Q 2 ' 



(ii) 



a s (Q 2 )/a s (Q 2 v ) if (Q 2 > Q 2 V ) 

1 if (Q 2 < Q 2 y) , 



where Q|=3.22GeV 2 , Qy=1.50GeV 2 , and the values of these special characteristic quan- 
tities are fixed by fitting data. 

Supposing the quark-hadron duality, the resultant sum rule for the mass of the doubly 
heavy baryon reads 



Ri 
R 



(12) 



with 



Ri 



1 



' T J (m d +m q ) 

T cond, 5/ 2 



ds P r\s)e- s/M * + B[nr d ' \q 2 )} + B[IT 



cond, 4/ 2 



(Q 2 )] 



+B[n: 

d 



)]+B[n cond, 6(g2)] 



dM 



2-R0 • 



(13) 
(14) 



B 



Here,m g (q=u, d, or s) denotes the masses of the light quarks, is the mass of the 
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diquark, Mb is the Borel parameter and so is the threshold cutoff introduced to remove 
the contribution of the higher excited and continuum states [30J. 



The perturbative contribution pi(s) and non-perturbative contributions 
B[Il™ nd ' dim (g 2 )] for Sqq/ and Qqq' in Eq. ffTB"]) are shown as follows: 

Pi{s) 



3 {m~ d - m- - s 



m, — m; 



+ s) 2 — Asm 



167TS 2 



B[n c r d ' 3 (g 2 )] = (M) 



m q (3M B + mjm 2 ^ m 



12M% 



q 'e m b , 



(15) 
(16) 



B[nr d ' 4,C (g 2 )] = (« S G 2 ) / dxe m b 



B[n c 1 ond,4,D (g 2 )] 



B[n cond, 4,E^ 2) ] 



lib j I'b r, 

a_j £. , r> 

x ( Q™2 



3x 2 - lOx + 10 



+ 



192ttM 2 192tt(x - l)x 2 M^ 



x 2 (-3a; 2 + lOx - 10) m 2 + (x - l) 2 (3s 2 - 7x + 3) 
1 



2x 2 (-3x 3 + 13x 2 - 21x + 11) 



384tt(x- 1) 2 x 3 M b 
+x\3x - 10)m^ + (x - l) 3 (3x 2 - 7x + 8m^ 



m 2 d m\ 



+ 



((x 2 — 4x + 3) m 2 — x 2 m d ) ((s — 1) 

1 



384tt(x - l) 3 x 4 M| 

2_2 J2„2\ 2 



'm- - arm 2 ,) 



(17a) 



9 2 
d_\ Q_ 

1 — x x 



{a s G 2 ) / dxe m b 
Jo 

-6 



+k„ + 2) 



512tt(x - 1) 2 3;M; 

/c„ ((x — l) 2 m 2 — x 2 m^) 



5127r(a; - 1)M 2 
x 2 m 2 d (—2xk v + k v — 2) + (x — l) 2 m 2 (2xK^ 

2 



(a,G 2 ) / dxe m b 
Jo 

-2) + Ak 2 v + 6^ + 4) 



512vr(x - l) 3 a; 2 M^ 
1 



9 2 
d_\ Q_ 

l — x X 



+ 



384tt(x - 1) 2 M 2 
1 



(17b) 

x{x 2 (6k v + 3) + x(Ak 2 — Ak v 
(x — l) 2 m 2 (x(6n v + 3) 



384tt(x - 1) 3 M B 
+4k 2 + 2k v + 1) - xm 2 d (x 2 (6k v + 3) + x {Ak 2 v - Ak v - 2) + 2k v ) 
1 



256n(x - l) 3 xM B 



(2/c„ + 1) (-ac (2x 2 - 3x + l) rn^m 2 , + x 3 m^ 



+ (x - l) 3 mj) 



(17c) 



B [nr d -v>] = -^w M - ( 3Mg(2 K „ +1 )- 8n »a 



»3 



g [n cond, 6(g2)] = g^giy (9Ml (2k v + 1) + 17m2) , (19) 

where the superscripts (C, D and E) on the left hand of Eq. (fI71) correspond to the labels 
in FigHJ 

For concision of the text, the detailed expressions of Rq are given in the Appendix. 



3 Numerical Analysis 

The numerical parameters used in this work are taken as 

(M) = -(0.254 ± 0.015GeV) 3 , a s (G 2 } = 0.07 ± 0.02GeV 4 , 
g s ($Ta ■ GV) = mjj(tftf>, a s (m) 2 = (2.1 ± 0.3) x 10~ 4 GeV 6 

m\ = 0.8 ± 0.2GeV 2 , m s = 0.14GeV, 



(20) 



m u ~m d = 0.005GeV, m c = 1.27GeV, 

m b = 4.19GeV, m Vc = 2.980GeV, 

mj/i , = 3.097GeV, m T = 9.460GeV. 

In our numerical analysis, we find that the effect of the two-gluon condensate is tiny, 
i.e., if we shift the contributions of two-gluon condensate a little, the final mass of the 
corresponding baryon hardly changes. We choose it as in Ref[28j : k v = 1.39. 

In practice, large uncertainty remains in the evaluation of baryon spectrum due to 
the input constituent quark masses. The mass of the baryon can be decomposed as 
mo + mq + AE, where tud, m 9 and AE are the masses of the heavy diquark, light-quark 
and the binding energy respectively. The binding energy is calculable within certain 
theoretical framework, whereas the quark, including the diquark, masses are usually not 
definite input parameters. We have to choose a reasonable strategy to determine the 
diquark masses, which influence more on the heavy baryon spectrum than light quark 
masses. It goes as follows. 
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There are few data on the doubly heavy baryons available, thus we cannot directly 
extract all necessary information from experimental data so far. Fortunately, the mass of 
S cc has been measured. Because of a simple symmetry argument, we know that the cc 
must constitute a spin-1 state, thus the cc diquark in S cc is an axial vector state. With 
the measured S cc mass as input, we determine ms cc y = 2.77 GeV. Then we still need to fix 
the masses of be and bb diqurks. Generally the effective potential between the two heavy 
quarks Q, Q' includes the Coulomb and linear confinement pieces. From the textbooks 
of the quantum mechanics, we know the solutions of the Schrodinger equations with the 
Coulomb potential proportional to ^ (the solution is the Lauerre polynomial) or the linear 
potential proportional to kt (the solution is the Airy function). The binding energy in 
the case where only the Coulomb piece exists is proportional to the reduced mass m of 
the QQ' system, instead, for the case where only the linear potential exists, the binding 
energy is proportional to (^) 1//3 if k is independent of m. Thus for the Schrodinger 
equation whose potential includes both the Coulomb and linear confinement pieces, the 
contributions from the two pieces compete and the dependence of the binding energy on 
the reduced mass is uncertain. Then we would like to invoke the data. 

As discussed above, m^ cc y and AE{ cc y can be directly extracted from the data, then we 
will use those values for the cc diquark and some proposed rules to fix m^ bc y, AE{ bc y, m^ bb y 
and AE{ bb y. Now let us make a plausible comparison of the quantities about diquarks 
with the corresponding mesons of the same flavor and spin structure. We can have AE^} 
and AE{ b iy from the relations: m^ b£ y = m b + m c + AE{ b£ y and m^ bb y = 2m b + AE^ bb y. 
Thus AE{ cS y and AE^ bb y are obtained as Mj/^ — 2m c and M T (is) — 2m b . 

The diquarks cc or bb are color-anti-triplet axial vector states, instead the meson J / ip 
or T(1S) are color-singlet vector states of cc and bb. The effective potentials between QQ 
and QQ only differ by a color factor, so that we may have 

AE{ cc y : AE{ c5 y = AE{ bc y : AE{ bE y = AE{ bb y : AE{ bb y, 

where the dependence of the binding energies on color and constituent masses may cancel. 
Then we obtain the binding energy for the axial vector bb. Since so far there are no data 
on B* available yet, we cannot determine AE{ bs y in the above scheme, but need to invoke 
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another way. Since be diquark is composed of c and b quarks, it is natural to think that 
an interpolation between cc and bb diquarks would be a good approximation for the be 
diquark, thus we write 

AE {bc} = -[AE m + AE {CC} ). 

Unlike the bb and cc diquarks, be diquark can be either an axial vector or a scalar. 
Now let us determine the mass of the scalar be diquark. The mass difference between 
spin-1 and spin-0 two-quark systems is due to the spin-spin interactions. Such interaction 
is proportional to l/(mQrriQi). Since there lack enough data for be mesons, let us first 
start with the charmmonia which are well measured and then generalize to the be mesons. 

The difference of the binding energies of J/ip and r? c is due to the spin-spin interaction 
between c and c, and besides a color factor related to the SU(3) Casimir factor which is 
| for a color singlet and | for a color- anti-triplet, the case for the cc diquark is the same. 
Thus we may write 

AE {cS} - AE [cS] = {M J/i} -M Vc ). 

Then using 

A-E{ c g} : AE[ c5 ] = AE{ bc y : AE[ bc ], 

we fix AE[ bc -\. 

With above analysis, the diqaurk masses which will be adopted in the following nu- 
merical computations are displayed as following: 

m J£ = 2.77GeV, m [cb] = 5.73GeV, m {cb} = 5.80GeV, m {bb} = 8.83GeV. (21) 

For choosing the proper threshold Sq and the Borel parameter Mj, there are two 
criteria. First, the perturbative contribution should be larger than the contributions 
from all kinds of condensates, and another is that the pole contribution should be larger 
than the continuum contribution [HI |29j [33]. On the other hand, the dependence of 
the evaluated masses of the doubly heavy baryons is rather unsensitive to variations of 
the Borel parameter in the Borel windows. For each baryon fortunately we can find 
an optimal Borel window where the two aforementioned criteria are satisfied and the 
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Table 2: The mass spectra of doubly heavy baryons. The "pole" stands for the contribu- 
tion from the pole term to the spectral density. The "cond" stands for the contribution 
from the condensate terms in the operator product expansion, where the threshold pa- 
rameter so takes its central value. Am is the energy gap between masses of other species 
of baryons and Hqq/ with the same diquark flavors. 



Baryon 


quark 


J P {J P D D ) 
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6.69 ± 0.08 + 5 2 cb 
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n* 
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°-' '-0.04 + °c6 
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{bb}s 
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1% 


0* 

"66 


{bb}s 


f(l + ) 


9.93l°;°I + 566 


0.13 


10.5-12.0 


(67-74)% 


2% 



results are almost independent of the Borel parameter after all. By the windows we 
obtain the masses of doubly heavy baryons. The dependence are shown in Figs. (J3H9]) , 
respectively. The numerical results are collected in Table [2] for various quantum numbers. 
For a comparison with other theoretical estimates on the baryon masses given in the 
literature, we also show those results in Table El The error bars are estimated by varying 
the Borel parameters, s° and the uncertainties in the condensates as well. It is noted 
that the uncertainty caused by introducing the diquark configuration is included in the 
diquark form factors (Eqs(10) and (11)). 

Note that inside the Tables (121131 . the mass of the baryon S cc with superscript "fit" is 
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Figure 3: Dependence of S cc and Q cc masses on the Borel parameter M|. The continuum 
thresholds s are taken as 3.9 2 , 4.0 2 , 4.1 2 GeV 2 for H cc , and 4.0 2 , 4.1 2 , 4.2 2 GeV 2 for ft cc , 
from down to up, respectively. We deliberately put two vertical lines denoting the chosen 
Borel window. 
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Figure 4: Dependence of S* c and Q* c masses on Borel parameter M|. The continuum 
thresholds s are taken as 4.0 2 , 4.1 2 , 4.2 2 GeV 2 and 4.1 2 , 4.2 2 , 4.3 2 GeV 2 , from down to up, 
respectively. We deliberately put two vertical lines denoting the chosen Borel window. 
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Figure 5: Dependence of "E ! cb and Q' cb masses on Borel parameter M\. The continuum 
thresholds s are taken as 7.0 2 , 7.1 2 , 7.2 2 GeV 2 and 7.1 2 , 7.2 2 , 7.3 2 GeV 2 , from down to up, 
respectively We deliberately put two vertical lines denoting the chosen Borel window. 
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Figure 6: Dependence of E cb and Q cb masses on Borel parameter M|. The continuum 
thresholds s are taken as 7.0 2 , 7.1 2 , 7.2 2 GeV 2 and 7.1 2 , 7.2 2 , 7.3 2 GeV 2 , from down to up, 
respectively. We deliberately put two vertical lines denoting the chosen Borel window. 
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Figure 7: Dependence of E* cb and masses on Borel parameter M|. The continuum 
thresholds s are taken as 7.1 2 , 7.2 2 , 7.3 2 GeV 2 and 7.2 2 , 7.3 2 , 7.4 2 GeV 2 , from down to up, 
respectively. We deliberately put two vertical lines denoting the chosen Borel window. 
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Figure 8: Dependence of E bb and Vt bb masses on Borel parameter M\. The continuum 
thresholds s are taken as 10. 7 2 , 10.8 2 , 10.9 2 GeV 2 and 10.9 2 , 11.0 2 , ll.l 2 GeV 2 , from down 
to up,respectively.We deliberately put two vertical lines denoting the chosen Borel window. 
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Figure 9: Dependence of and Ql b masses on Borel parameter M|. The continuum 
thresholds s are taken as 10.4 2 , 10.5 2 , 10.6 2 GeV 2 and 10.5 2 , 10.6 2 , 10.7 2 GeV 2 , from down 
to up, respectively. We deliberately put two vertical lines denoting the chosen Borel 
window. 




Table 3: Comparison with other theoretical resutls and the experimental data (if avail- 
able). The quantities are in GeV. 
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J P (Jn°) 


Our work 
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Exp. [9] 
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2 v x ! 
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2 V- 1 - ) 
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{bb}s 
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2 V ± / 
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taken as inputs to obtain the mass of diquark m{ cc } and then the masses of other baryon 
states in the table are predicted. 

As indicated in above the choice of the diquark masses is based on our postulate 
about the binding energy, so this strategy would certainly bring up some theoretical 
uncertainties. To explicitly show how the diquark mass influences the spectrum of doubly 
heavy baryons, let us shift the corresponding diquark masses by 0.1 GeV, and we find that 
the uncertainty of the baryon mass lies within 0.047 ~ 0.064 GeV. One can be convinced 
that the uncertainty should be no more than 10% as we change the diquark mass within 
a reasonable range. Moreover, in the forth column of Tabled we put a term Sqq' following 
the predicted mass to manifest a possible error. In next column of this table, we list the 
gaps (Am) among the concerned baryon masses where the uncertainties cancel out, and 
hence may make more senses. That means the predictions in this work, especially on the 
mass gaps, are experimentally testable. 

4 Summary and Conclusions 

In this work, the masses of various doubly heavy baryons have been studied in terms 
of the QCD sum rules where the diquark structures are priori assumed. In the calculation 
we keep the contributions of the condensates up to dimension six in OPE. Our results, 
in certain tolerance, are in accordance with the theoretical predictions via other models. 
Especially, it is worth pointing out that our results are reasonably consistent with that 
calculated in the QCD sum rules without assuming diquark structures. 

In the calculation, an effective coupling between diquark and gluon which was phe- 
nomenologically introduced is adopted. The form factor at the effective vertex indeed 
manifests an inner structure of the diquark. But as the diquark is viewed as an indepen- 
dent degree of freedom, this factor performs as an ad-hoc parameter in the given theory 
and it plays a role just as the quark or gluon condensates in the QCD sum rules which 
were obtained either from an underlying theory (such as the value of the gluon conden- 
sate could be obtained from the dilute gas approximation of instantons) or by fitting data 
(such as the value of the quark condensate might be gained by fitting the pion decay 
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constant). 



Our results imply that the structure of a heavy diquark and a light quark is indeed 
a reasonable configuration for the doubly heavy baryons. The Large Hadron Collider 
(LHC) which has already begun running, even at lower energy (7 TeV) and luminosity , 
will provide a large database of doubly heavy baryons. Once enough data are available, 
one can further analyze the doubly heavy baryons of various flavors and spins. Comparing 
our theoretical predictions on their mass spectra with the data, will not only enrich our 
knowledge on the underlying theory, i.e. the low energy QCD, but also further investigate 
the diquark structure and applicability for dealing with the processes such as production 
and decay of the doubly heavy baryons. 
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Appendix 

The perturbative contribution p2(s) and nonperturbative contributions B[il2° nd ' dim (? 2 )] 
for Sqq/ and flQQ> in Eq.f fT3"j) are shown as follows: 
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The perturbative contributions Pi(s) and nonperturbative contributions B[Il^ ond ' (q 2 )] 
for Hgg, and Qqqi in Eg. ({TBI are shown as follows: 
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The perturbative contributions pi{s) and nonperturbative contributions B[H^ ond ' dim (g 2 )] 
for E'qq, and Q'qqi in Eq. ffTBl are shown as follows: 
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